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Abstract 

The question of conductivity is revisited. Using the total mo- 
mentum shift operator to construct the perturbed many-body Hamil- 
tonian and ground state wave function the second derivative of the 
ground state energy with respect to the perturbing field is expressed 
in terms of the one and two-body momentum densities. The distinc- 
tion between the adiabatic and envelope function derivatives, hence 
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that between the Drude and superfluid weights can be introduced in 
a straightforward manner. It is shown that a discontinuity in the 
momentum density leads to a contribution to the Drude weight, but 
not the superfluid weight, however a (5-function contribution in the 
two-body momentum density (such as in the BCS wave-funtion) con- 
tributes to both quantities. The connection between the discontinuity 
in the momentum density and localization is also demonstrated. 

To distinguish between conductors and insulators an expression for the 
frequency-dependent conductivity was derived by Kohn [1]. The DC conduc- 
tivity (Drude weight) corresponds to the strength of the (5-function peak of 
the conductivity at zero frequency. The Drude weight is often expressed [H [2] 
in terms of the second derivative of the ground state energy with respect to a 
phase associated with the perturbing field. This phase has the effect of shift- 
ing the momenta of the system. Scalapino, White, and Zhang (SWZ) [3l H] 
have pointed out that taking the derivative with respect to the phase is 
ambiguous: if the derivative is defined via adiabatically shifting the state 
which is the ground state at zero field, then the Drude weight results. In the 
presence of level crossings the adiabatically shifted state may be an excited 
state for finite perturbation. The superfluid weight is obtained if the deriva- 
tive corresponds to the "envelope function", i.e. the ground state for any 
value of the perturbation. SWZ also state that nonadiabatic crossings occur 
infinitesimally close to zero field if the dimensionality is greater than one. 

In this paper this question is revisited. Based on the total momentum 
shift operator [5] the perturbed Hamiltonian and ground state wavefunction 
are explicitly constructed. This operator plays an important role in con- 
structing the total position operator for many-body systems [SI El C] • The 
second derivative of the ground state energy with respect to the perturbing 
field is then expressed in terms of the one and two-body momentum den- 
sities. It is then shown that the adiabatic and envelope derivatives can be 
distinguished by varying the length scale associated with the total momen- 
tum shift operator, which is also the length scale of the perturbing field. 
When this length scale is assumed to be the same as the size of the system 
then ^ ^2^-* is proportional to the superfluid weight, if this length scale is 

assumed to be much larger than the system size than ^ ^2^'' corresponds to 
the Drude weight. For continuous one and two-body momentum densities 
both quantities are zero. If the one-body momentum density is discontinu- 
ous then the Drude weight is finite, but the superfluid weight is zero, and if 
the two-body momentum displays a 5-peak (Cooper pairing) then both the 
Drude and superfluid weights are finite. Hence insulators, metals, and super- 
conductors can be distinguished. While a discontinuous momentum density 



2 



being a sign of conduction is a well-known result of many-body theory [S] 
and plays an important role in the Landau theory of Fermi liquids |S1 E], 
the foundations of the latter are distinct from those for the conductivity put 
forth by Kohn pLj. In this work the finiteness of the Drude weight and the 
discontinuity in the momentum density are shown to coincide. Moreover, 
it is also demonstrated that the localization tenet suggested by Kohn pQ, 
namely that a system localized(delocalized) in the many-body configuration 
space is insulating(metallic), is also equivalent to the absence(presence) of a 
discontinuity in the momentum density. Hence the Landau theory of Fermi 
liquids and the localization theory of Kohn are placed on the same theoretical 
footing. 

We consider a system of interacting fermions whose Hamiltonian is peri- 
odic in L. We will assume that the ground state is also periodic in L (i.e. 
$ = 0). This leads to no loss of generality, since if the ground state is at a 
finite the Hamiltonian can be shifted. We wish to write the Hamiltonian 
for such a system. We first write 

^ = ^(M^)};{£?}) (1) 

where g{k) are continuous functions of k and c^^"* denote creation and anni- 
hilation operators of particles at wave-vector k. This Hamiltonian includes 
only states which are periodic in L. Due to the periodicity the spacing of 
the points on which the momenta are represented is A/c = 27r/L. H is not 
the full Hamiltonian of the system, since the states with twisted boundary 
conditions (which correspond to A;- vectors which fall between the grid-points) 
do not appear as eigenstates. To include them we write 

n{a)=H{{g{k + a)y,{cf^^}). (2) 

Here all the k vectors have been shifted by a, however, the spacing of the 
fc-vectors is unchanged. The full Hamiltonian can be written 

'H^ = -L r da-H(a). (3) 

This Hamiltonian is the full Hamiltonian in the sense that the system itself is 
periodic in L, however states of all boundary twists are included. Ht is block- 
diagonal, since Hamiltonians with different values of a correspond to different 
Hilbert spaces. In the limit L oo Ht becomes the full Hamiltonian of the 
infinite system. 

To stress this point one can consider the Hubbard model for a system 
with size L with Hamiltonian written in reciprocal space, 

Hnub = ^kanka + UY1 cl^cl,^Ck+qfCk' -gl- (4) 
kcr kk'q 
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The eigenstates of this Hamiltonian are periodic in L. The spacing of the 
fc- vectors is Ak = 271 /L. The shifted Hubbard Hamiltonian 

(5) 

ka kk'q 

has eigenstates with twisted boundary conditions, however, the Hamiltonian 
still corresponds to a system periodic in L, as the spacing between the k- 
vectors is still A A; = 2tt/L. 

It is expedient to introduce the total momentum shift operator 

where X = the sum of the positions of all the particles, and which 

has the property that [5] 

We extend U{2tt/L) to lengths nL with n integer. Then momentum shifts 
to states with twisted boundary conditions on L are also included. Taking 
the limit — t- oo we can write 

f/(7)?i(«)f/(-7) = nigik + a)}; {cIVt))' («) 
for arbitrary 7 thus 

U{^)nTUi-l) = ^ £daH({(7(fc + a + 7)};{cgj). (9) 

The transformed Hamiltonian defined in Eq. (|9]) has the same eigensystem 
as Ht- The transformation merely shifts the block diagonal Hamiltonians 
which comprise Ht- 

The linear response of a system with periodic boundary conditions can be 
cast using the total momentum shift. We assume that the system of interest 
has a Hamiltonian of the form 

H = J2^knk + Hi, (10) 

k 

where Hi denotes an interaction diagonal in the coordinate representation. 
This Hamiltonian includes the ground state, which is also periodic in L. For 
the ground state wavefunction we assume the form, 

\^m= y: mi-ikN)ci...ci\o), (n) 

fcl,...,fcjv 
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which is the most general for fixed particle number. 

The usual way to introduce a static vector potential Ax is to multiply 
the hopping parameters with a phase factor. In this case the k vectors are 
shifted as — )■ + $ with $ = A/ he, leading to 

Hi<^) = J2ek+^nk + Hi. (12) 

k 

To arrive at Eq. (|T2|) one can also use the total momentum shift opera- 
tor on the total Hamiltonian constructed from H, and shift indices as was 
done to obtain Eq. In the same way one can obtain the wavefunction 
corresponding to the shifted H{^), 

|v^($))= 5: ^(fci + $,...,fc;v + $)4...4|0). (13) 

fci,...,fcjv 

The criterion for the DC conductivity and the superfiuid weight can both 
be written [21 El H] in the form 

D = (14) 

While the Drude weight and the superfiuid weight quantities correspond to 
different perturbations, the expression for these quantities coincides, since 
in the above expression $ = 0, hence the explicit dependence on the vector 
potential, which gives rise to the distinction, is neglected. Taking advantage 
of the Hellmann-Feynman theorem D can be expressed as 

1 i,^, ,,92iJ(0),^, ,9^(0),9/7(0),^, ,,9i7(0),9^(0)/ 

(15) 

The reason that both the Drude and superfiuid weights can be written in 
this form is due to the fact that Eq. (|T5l) is a linear response expression in 
which the effect of the perturbing field is set to zero. 

The derivatives with respect to $ of the Hamiltonian can be made to 
correspond with derivatives with respect to the momenta, i.e. it holds that, 

and 

-9^ = E^^n,, (17) 
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The derivative with respect to k is ambiguous [21 II]- For a finite system with 
size L the summation in Eqs. ( 1T6|) . ( TT71) . and (fT8|) is defined on grid points 
separated by 27r/L in reciprocal space. Thus one way to define the derivatives 
is using these grid points (for example the finite element definition). 

The total momentum shift extended to length nL extends the Hilbert 
space, hence the derivatives can also be defined using the extended states on 
the finer grid 27c /{nL). Note that the summations in Eqs. f|T6|) . f|T7|) . and 
f|T8|) are still defined on the grid 211 /L. When the thermodynamic limit is 
taken is a continuous function, hence this distinction between grids causes 
no ambiguity in the application of Eqs. ffTB]) and f lT7|) . The wavefunction, 
however, can be discontinuous, and, as discussed below, this leads to conse- 
quences. Using Eqs. (|T6|) . ( IT7|) . and (|T8i) one can show that 



where Uk and n\ y denote the one and two-body momentum densities in the 
ground state, defined as 

^fc = E E l^(fci,-,MI' (20) 

i fci,...,fcjY 

and 

4'1' = E E mu-.k^)?. (21) 

Eq. (fT9|) is arrived at by using Eqs. (|T6|) . ( 1T7|) . and (fTSi) . and the identity 

{0\ck^...Ck,hkcl^...cljO) = E'^fc.fc- (22) 

i 

For the case n = 1, we replace the derivative in Eq. (HM by 

duk nk+2n/L - nk , , 

dk ^ 27r/L • ^^^^ 

This definition corresponds to the "envelope function" definition of SWZ |21 
H]. To see this consider the system at $ = and $ = 2it/L. The ground 
state at $ = is of the form in Eq. ([TT]), at $ = 2tt/L it is Eq. ([13]), no 
longer the ground state in general. For both $ = and $ = 2%/ L the ground 
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state density is given by ra^. In Eq. fl2^ the function (corresponding to 
the ground state) is used in both cases. When the thermodynamic hmit 
{L —7- oo) is taken the first two terms in Eq. (|T9i) cancel due to partial 
integration resulting in 

i^("=^) = A / dkdk'^^. (24) 



TT^ J dk dk' 

This quantity integrates to zero, due to the periodicity of the Brillouin zone, 
unless, as discussed below, pairing occurs in the two-body density. These 
arguments allow association of Z}("=^) with the superfluid weight. 

We now consider the implications of the different properties of the deriva- 
fives for n = 1 and n — )■ oo. For segments for which nk and ^, are con- 
tinuous the two definitions of the derivatives (based on the spacing 27r/L 
vs. 27r/(nL)) coincide, however this is not true when either densities are 
discontinuous in k. While on the larger grid 27r/L a discontinuity in these 
quantities leads to a divergence, on the grid 27r/(nL) the discontinuity does 
not occur when the derivative at the fc-grid points is evaluated and the limit 
n — )• oo is taken first, and the derivative is defined as adiabatically shifted. 

As an example one can consider a Fermi sea, for which the term depend- 
ing on the two-body density does not contribute since there are no corre- 
lations between momenta. When a phase is applied the energy levels and 
the momentum densities are shifted as — )■ ek+^, — )■ nk+^. If the phase 
$ ^ 2tt/L, and the ground state of the new Hamiltonian is used in defining 
the derivative ("envelope function"), then the discontinuity contributes to 
the derivative, since if is the last filled state near the discontinuity, then 
nk+2n/L will be the first unfilled one. However, for small $ (which corre- 
sponds to the limit n — t- oo) if corresponds to the last filled state then 
nA;+$ does not change. Excluding the discontinuities (which are relevant to 
the second term in Eq. (fT9|) ) from the partial integral leads to 

= ^An,,^, (25) 
In ok 

where the discontinuities are assumed to he at k = ±kF (Fermi wave vector). 
When spin is included then each spin component will contribute a term of 
the form in Eq. fl25p . For this reason we associate the quantity Z)*^"~^°°) with 
the Drude weight. 

To explore the connection between conduction and the discontinuity in 
the momentum density further we consider the quantity 



n{y)= {mivm 



J2 ^* iki + y,..;kN + y)ip{ki,...,kN] 

i,...,A;jv 



(26) 
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The quantity (L^/ (27r^)Re lnn(27r/L) was suggested by Resta and Sorella as 
a criterion of localization. As a result of Kohn's hypothesis [T] localization 
is also a criterion to distinguish conductors from insulators. If the wavefunc- 
tion ip{ki, ...^k^) is a continuous functions of its arguments then n(2'7r/L) 
approaches unity in the limit of large system size. The functions Uk and n(?/) 
are then continuous, corresponding to insulation. When is discontinuous 
then the magnitude of the wavefunction ip{ki, k]\f) is also discontinuous. 
In the following we assume that the magnitude of tplki, k^) is discon- 
tinuous but its phase is not. Since ip{ki, k^) describes indistinguishable 
particles, the discontinuity has to occur as a function of any of its arguments. 
Moreover, on physical grounds we anticipate that this discontinuity occurs 
at the Fermi wave-vector. The effect of the discontinuity can be assessed by 
considering the difference 

n(0)-n(e) = \j^Jdki.dkNr{ki,...,kN)^{ki,...,kN) 

- I dki.dkN^*{ki + e,...,kN + e)tp{ki,...,kN)\. (27) 

where e denotes an infinitesimal and the thermodynamic limit was taken. The 
integrands in the first term and the second term will cancel for regions where 
the coefficient ipi^ki, ki\i) is continuous. The contribution of a discontinuity 
Sit kp will be of the form 

p{kp+; kF+) + pikp-; kp-) - p{kF+; kp-) - pikp-] kp+), (28) 

where p{k; k') denotes the one-body density matrix in fc-space. Rewriting in 
a natural orbital representation this contribution takes the form 

E^i[7r(^F+)7i(fcF+)+7;(^F-)7*(fcF-)-7;(fcF+)7*(fcF-)-7;(^F-)7^(^F+)], 

(29) 

(with < qi < 1 and 7j(A;) denoting the natural orbitals) which under the 
assumption of a continuous phase is a positive quantity. Since this is also the 
case for 11(0) — n(—e) it follows that for discontinuous coefficient ip{ki, k^) 
the function Il{y) will contain a 5-function contribution at the origin. These 
results coincide exactly with the results of Resta and Sorella [7] where a 
function of the quantity |n(27r/L)| is suggested as a criterion of localization 
and conduction. 

To understand the effect of pairing we study the BCS wavefunction 

\^BCs) = Ui^k + Vkclcl,^)\0). (30) 

k 

We assume a BCS Hamiltonian with constant coupling between Cooper pairs. 
Calculating the properties of this wavefunction requires generalization to 
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include spin and variable particle number of Eq. f[T^ which presents no 
difficulty. Since the one-body density of the BCS wavefunction is continuous 
the ffist two terms cancel by partial integration when the thermodynamic 
limit is taken. Thus we are lead to consider the last term only, which depends 

(2) 

on the two-body momentum density n^^/. This quantity can be broken up 
into components with parallel and anti-parallel spins. The parallel spin two- 
body density is again continuous, hence does not contribute. The two-body 
density when the spins are anti-parallel gives 

- I fxk)f{k') k' ^ -k, ^^'^ 

with 

fik) = „ • (32) 



Explicit calculation for the BCS wavefunction then yields for the thermody- 
namic limit 

The first term arises since n^^a-k-a = ""-fco-; i-e- due to Cooper pairing. Due 
to the continuity of Uka the last term is zero. Partial integration then results 
in 

1 V- f ,.d-tt 



(33) 

Since the function f{k) is continuous this result holds for both n = 1 and 
n — )■ oo. The result that the second derivative of the "envelope" function 
of the ground state energy is finite for a superfiuid and zero for a normal 
metal was obtained for the case of a ring with finite thickness by Byers and 
Yang [10]. 

SWZ have also shown [1] that for dimensions higher than one the first 
non-adiabatic crossing occurs at zero field when the thermodynamic limit is 
taken. This leads to a distinction between evaluating ^ Q^f' first and then 
taking the thermodynamic limit or vice versa. In generalizing the formalism 
presented here to higher dimensions one has to consider that the differential 
operators in the superfiuid and Drude weights operate in one particular direc- 
tion (that of the perturbing field). If the thermodynamic limit is first taken 
in the direction perpendicular to the perturbing field, then the discontinuity 
can "disappear" . For example, a two-dimensional non-interacting system at 
half-filling has a discontinuous momentum density, Uk^^ky, but the function 
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f{kx) = I dkyUk^ ky is a continuous function. However, the definition of tlie 
derivative corresponding to the case n — t- oo resolves this ambiguity. In 
that case irrespective of the order of hmits the discontinuity will be excluded 
from the integration, as argued above for the Fermi sea. Moreover, as shown 
above, the discontinuity in the momentum density contains exactly the same 
information as the localization order parameter of Resta and Sorella [7], a 
quantity which is also insensitive to dimensionality. 

In conclusion the second derivative of the ground state energy with respect 
to a perturbing field (vector potential) at zero field was derived and shown 
to be an expectation value over the one and two-body momentum densities. 
A length scale associated with the perturbation was defined, and through 
it states with twisted boundary conditions were introduduced, allowing for 
the possibility of defining the adiabatic derivative (Drude weight) and the 
derivative of the ground state energy envelope function (superfiuid weight). 
The resulting expression for the Drude weight is not the zero frequency limit 
of an quantity based on time-dependent perturbation theory. The Drude 
weight is finite in the presence of discontinuities in the wavefunction (which 
correspond to discontinuities in the momentum densities), as well as due to 
BCS pairing. The superfiuid weight is not sensitive to discontinuities in the 
momentum densities, but is finite in the presence of BCS pairing. It was 
shown that a localization quantity suggested by Resta and Sorella [7] based 
on a tenet of Kohn [IJ contains the same information as the discontinuity 
in the momentum density. Thus the connection between the localization 
hypothesis of Kohn [T] and the criterion of metallicity in the Landau theory 
of Fermi liquids is established. 
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